IDEAS FOR DESIGN 



Circle S2Q 



Protect Your Data 
With Forward Error Correction 

S.K.SHENOY 

Naval Physical & Oceanographic Lab., Thrikkakara PO, Kochi 682 021, 
India; e-mail: kedar@satyam.net in 



Multiprocessor architectures 
have become commonplace in 
embedded applications these 
days. The processors often are inter- 
connected via simple asynchronous 
serial links. These links employ 
UARTs that don't offer any built-in 
error control other than simple one- 
bit error detection using parity. 
Many designers (except for commu- 
nication specialists) are unaware that 
it's not difficult to incorporate reason- 
ably powerful error correction in 
software. 

T he performance overheads are 
\ i tolerable in most applications. In 
fact, this error-correction method can 
be reduced to simple table lookups. 

The method described here, based 
on the well-known Hamming code 
technique, enables correction of sin- 
gle-bit errors. The method also can be 
extended to burst errors as explained 
below. This forward-error-correction 
method is useful in applications where 
data integrity is important. In addi- 
tion, it's beneficial in simplex systems 
or long transmission delay systems in 
which error detection, acknowledge- 
ments, and retransmission aren't fea- 
sible. One good example is an under- 
water acoustic data-communication 
system. 

The technique illustrated here uses 
4-bit data characters. Three parity 
bits are added to this data character 
to form the 7-bit character that's 
transmitted. This character size is 
very convenient, since it's supported 
by most UARTs and other communi- 
cation devices. 

Let the Hamming-encoded data 
3 be numbered as shown below, 
wnere Bl is the LSB: 

B1B2B3B4B5B6B7 
P1P2D1P3D2D3D4 

The corresponding bits P n and D n 
indicate how this character is assem- 



bled. PI to P3 are the inserted parity 
bits, which are computed as shown be- 
low, while Dl to D4 are the data bits. 
Note that the P n bits occur at posi- 
tions corresponding to powers of 2. 
Representing each B n in terms of 
B m s, which correspond to powers of 2, 
one can write: 

B3 = B1 + B2 
B5 = Bl + B4 
B6 = B2 + B4 
B7 = Bl + B2 + B4 

Bl occurs in the expansions for B3, B5, 
and B7; B2 in B3, B6, and B7; and B4 in 
B5, B6, and B7. 

Using this information, the Ham- 
ming coding requires that Bl be rep- 
resented as the parity of bits B3, B5, 
and B7, and so on. Thus, the parity 
bits are computed in terms of the data 
bits as: 

PI = Dl + D2 + D4 
P2 = Dl + D3 + D4 



P3 = D2 + D3 + D4 

where + indicates the Exclusive-OR 
operation. 

On the receive side, the values of 
the PI, P2, and P3 bits are recom- 
puted and compared with the received 
values of PI, P2, and P3. The result is 
then assigned to R1-R3 such that if 
the received and computed PI match, 
then Rl = 0. Otherwise, Rl = 1; like- 
wise for R2 and R3. If Rl-3 is treated 
as a 3-bit number, it will indicate the 
bit in error as shown (Table 1). 

The erroneous bit is then corrected 
by inverting it. The attached C pro- 
gram implements the technique de- 
scribed (see the listing, nextpage). Be- 
cause 8-bit data is typically used, the 
program takes each 8-bit character, 
splits it into two 4-bit nibbles, then en- 
codes and transmits each nibble as a 7- 
bit character. The 8-bit character is re- 
assembled on the receive side after 
error correction. 

This technique also can be ex- 
tended to correct burst errors. For 
correction of burst errors up to seven 
bits in length, seven characters are as- 
sembled and stacked as shown to form 
ablock (Table 2). 

The vertical slices corresponding to 
each bit position (Al to Gl) are then 
taken to form a 7-bit character and 
transmitted without any further en- 
coding. All seven such characters are 
thus transmitted. On the receive side, 



TABLE 1: BIT-ERROR INDICATION 
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B1 in error 
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B2 in error 
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B3 in error 
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B4 in error 
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B5 in error 
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B6 in error 
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B7 in error 








TABLE 2 


: BLOCK FORMATION 


Chan 


A1 


A2 


A3 


A4 


A5 


A6 A7 


Char 2 


B1 


B2 


B3 


B4 


B5 


B6 B7 


Char 3 


C1 


C2 


C3 


C4 


C5 


C6 C7 


Char 4 


D1 


D2 


D3 


D4 


D5 


D6 D7 


Char 5 


E1 


E2 


E3 


E4 


E5 


E6 E7 


Char 6 


F1 


F2 


F3 


F4 


F5 


F6 F7 


Char 7 


G1 


G2 


G3 


G4 


G5 


G6 G7 
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* Protect Your Data Using Forward Error Correction 

* S.K.Shenoy, NPOL, Kochi, India 
*/ 

#include <stdio.h> 
#include <stdlib.h> 
#include <conio.h> 
#include <dos.h> 

typedef unsigned char byte; 

r Returns 1 if ith bit of the character 

* 'value' is set, else returns 7 
byte BitSet( byte value, byte i ) 

{ 

byte bit = 0x1 ; 

bit = bit « i; 

if ( ( value & bit ) = bit ) 

return( 1 ); 

else 

return( ); 

} 

r Sets a specified bit of a character to 1 or 7 
void SetBit( byte *Val, byte BitPos, byte BitVal ) 



/•Reset bit 7 



{ 

if ( BitVal == ) 
*Val = *Val & -( 0x01 « BitPos ); 
else 

*Val = *Val I ( 0x01 « BitPos ); 

} 

r Splits the input 8-bit character into two nibbles, encodes each nibble 

into 7-bit hamming codes by adding the 3 parity bits. 7 
void GenHammingCodes( byte chr, byte *ch1 , byte "ch2, byte CodeTab[ ] ) 



{ 

*ch1 = ( chr & OxOF ); 
*ch2 = (chr&0x0F0)» 4; 
*ch1 = CodeTab[ *ch1 ]; 
*ch2 = CodeTab[ "ch2 ]; 
} 

r Inverts the specified bit 7 

void lnvertBit( byte *Val, byte BitPos ) 

{ ....... 

if (Bitset( *Val, BitPos )) 
SetBit( Val, BitPos, ); 
else 

^SetBit(Val, BitPos, 1); 

/* Inverts the specified bit 7 
void CorruptBit( byte *Val, byte BitNo ) 
{ 



r Get LS Nibble of data 7 
r Get MS Nibble of data 7 



I 

byte i; 



for(i = 0;i<16; i++) 

CodeTab[ i ] = ComputeHammingCode( i ); 

} 



I" Generates the 128-entry Decode Table of 7-bit c 

for "corrupted" codes 0000000 to 1 1 1 1 1 1 1 7 
void GenerateDecodeTable( byte DecodeTab[ ], byte CodeTabf ] ) 
{ 

byte i, DataNibble, ComputedCode, CorrectCode, CorruptBitlndex; 

/* i represents the 128 possible 'corrupt' 7-bit codes 7 
for(i = 0;i<128;i++) { 

DataNibble = ( i & 0x78 ) » 3; /* Get only the 4-bit data nibble 7 

f Get correct code for 'corrupt' code 7 

ComputedCode = CodeTab[ DataNibble ]; 

CorruptBitlndex = ( i & 0x7 ) A ( ComputedCode & 0x7 ); 

CorrectCode = i; /* Correct the 'corrupt' code 7 

if( CorruptBitlndex ) /* Non-Zero value = the bit number which is corrupt7 

{ 

f The following if and else if statements compensate for 
the exchange of the positions of P3 and DO bits in forming 
the character/ 

if( CorruptBitlndex = 3 ) CorruptBitlndex = 4; 

else if( CorruptBitlndex == 4 ) CorruptBitlndex = 3; 

lnvertBit( &CorrectCode, CorruptBitlndex-1 ); /* Correct the error 7 

} 

r Fill the table 7 



lnvertBit( Val, BitNo ); 

) 

/* Encodes the input 4-bit value into 7-bit hamming code 

by adding the 3 parity bits. 7 
byte ComputeHammingCode( byte chr ) 



DecodeTab[ i ] = CorrectCode; 
} 
} 



r Corrects the errors in the two receive-encoded characters, strips the 

parity bits and re-assembles the original 8-bit character. 7 
void ReAssembleChar( byte "ch1 , byte "ch2, byte "ch, byte DecodeTab[ ] ) 
{ 

byte c1,c2; j 
d = DecodeTab[ "ch1 ]; /* Correct the first code word 7 
c2 = DecodeTabj *ch2 ]; /* Correct the second code word 7 
r Combine the two data nibbles 7 
*ch = ((d »3)&0xF)l((c2« 1 )&0xF0); 

void main( void ) 
{ 

byte val ■ 'a'; 

byte *ch; 

byte vaM,val2; 

Irit BitNo; 

byte CodeTab[16]; 

byte DecodeTab[ 128]; 

r Generate the Encode and Decode tables 7 
GenerateCodeTable( CodeTab ); 
GenerateDecodeTable( DecodeTab, CodeTab ); 



byte code, pO, p1,p2; 
code = ( chr &0xF)« 3; 



r Get LS data Nibble in 4 MS bits 
of 7-bit Code word 7 



r Forever Loop 7 
T Key pressed 7 
I" Read user data 7 

r Exit if Escape key pressed 7 



r Compute the three parity bits 7 

pO = ( BitSetf code, 3 ) A BrtSet( code, 4 ) » BitSet( code, 6 )) & 0x1; 
p1 = ( BitSet( code, 3 ) A BitSet( code, 5 ) A BitSet( code, 6 )) & 0x1; 
p2 = ( BitSet( code, 4 ) A BitSet( code, 5 ) A BltSet( code, 6 )) & 0x1 ; 
r Insert the parity bits in the 3 LS bits of code word 7 
SetBit( &code, 0,p0); 
SetBit( &code, 1, pi ); 
SetBit( &code. 2, p2 ); 
retum( code ); 



} 



T Generate 16-entry Coding Table of 7-bit 

for data nibbles 0000 to 1111 7 
void GenerateCodeTable( byte CodeTab[ ] ) 



•A 



C3 



while ( 1 ) 
{ 

H ( kbhit( ) ) 
{ 

val = getch( ); 
If (val = 27) 

exit( ); 

}> , j • • ' 

r Generate the 2 Hamming codes 7 
GenHammingCodes( val, &val1, &val2, CodeTab ); 
prlntf( 'val = %x vail = %x val2 = %x\n", val, vail , val2 ); 

printf( "Enter bH position to be corrupted ( -> 6 ):" ); 
fscanf( stdin, "%d", SBitNo ); 

CorruptBlt( &val1 , (byte) BitNo ); r Corrupt the specified bit 7 

CorruptBit( &val2. (byte) BitNo ); 

printf( "val = %x vail = %x val2 = %x\n", val, vail , val2 ); 

f Correct & re-assemble char 7 

ReAssembleChar( &val1, &val2, ch, DecodeTab ); 

printf( "ch = %x\n', *ch ); 

} 

> 1 8 «H» i • q .. , 



j. nc cii ux correctr 
carried out on each reassembled char- 
acter (Char 1-7). 

In the event of burst errors, at 
most one bit of each character (Char 1 
to Char 7) will be corrupted and can 
be corrected. For this format, data 
must be available in blocks of seven 
characters, which may not be the case 
in many applications. In such in- 
stances, the first character, Char 1, 
can be used to indicate the actual 
number of data characters that follow 
in the block. The rest of the charac- 
ters, which are dummy characters in- 
serted to make up the seven-charac- 
ter block, can be ignored. The block 
data count also is protected by the er- 
ror correction. Note that there are a 
number of ways of marking packet 
boundaries in asynchronous commu- 
nication systems. 

The demo C program takes the data 
character input by the user and splits 
it into two 4-bit nibbles. Each nibble 
then is encoded by adding three parity 



hamming code. A single-bit 
error can be simulated on these codes 
at a user-specified bit position by in- 
verting the bit concerned. On the re- 
ceive side, the error is corrected and 
the original 8-bit character is reassem- 
bled from the two data nibbles ex- 
tracted from the encoded characters. 

The encoding table is generated by 
applying the encoding algorithm to 
each of the 16 4-bit patterns (0000- 
1111), resulting in a 16-entry table 
with each entry containing a 7-bit en- 
coded value. The encoded value for the 
data nibble xxxx is given by the 7-bit 
pattern at the xxxx table position. 

The decoding table is generated by 
applying the decoding algorithm to 
each of the 7-bit patterns (0000000- 
1111111), representing the possible 
received values of encoded data. This 
creates a 128-entry table with 7-bit en- 
try values. If the received bit pattern 
is xxxxxxx, the 7-bit entry at index 
xxxxxxx represents the corrected 
data character. 



and then used for all subsequent en- 
coding/decoding activity. 

The program was compiled and 
tested in Borland C 4.5 on a PC. The 
encode and decode routines can be 
easily ported to almost any C environ- 
ment, including cross compilers for 
microcontrollers. The encoding and 
decoding tables can be generated and 
stored in PROM memory. 

A note on the program listing: The 
data and parity bits have been inter- 
leaved in the above discussion for sim- 
plicity. This isn't a requirement. In the 
program listing, the four data bits and 
three parity bits are conveniently 
packed together as shown below to 
form the 7-bit data: 

B1B2B3 B4 B5 B6B7 
P1P2P3D1D2D3D4 

The only effect of this rearrange- 
ment is that the "Bit In Error" desig- 
nations shown in Table 1 must be ad- 
justed to indicate the new bit positions. 
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T 



he circuit presented can be helpful 
in applications where the desired 
output is a dc offset voltage of an 



ac signal. Such an application would be 
the design of an electronic watt/watt- 
hour meter. For example, when calcu- 



lating active ac power, Vp*sin(cot) 
(proportional to voltage across the 
load) and V c *sin(cot + \|0 (propor- 
tional to the current through the 
load) are multiplied, resulting in 
K*V P *V c *cos(v) [dc offset] and 
K*V P * V c *cos( 2cot + i|f) [ac harmonic]. 

Because the required output signal 
is just the dc term, the ac term must 
be eliminated. This type of function 
can't be performed by conventional 
low-pass filters. However, the simple 
and inexpensive circuit shown accom- 
plishes the job (Fig. 1). 




1 . This simple op-amp circuit is used to extract the dc offset voltage from an input signal containing low-frequency ac powerline components. 



